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ABSTRACT 

This paper presents a new algorithm, TRIMOR, to analyse multi-order spec- 
tra of triple systems. The algorithm is an extension of TRICOR, the three- 
dimensional correlation technique that derives the radial velocities of triple 
stellar systems from single-order spectra. The combined correlation derived 
from many orders enables the detection and the measurement of radial veloc- 
ities of faint tertiary companions. The paper applied TRIMOR to the already 
available spectra of HD 188753, a well known triple system, yielding the radial 
velocities of the faintest star in the system. This rendered the close pair of 
the triple system a double-lined spectroscopic binary, which led to a precise 
mass-ratio and an estimate of its inclination. The close-pair inclination is very 
close to the inclination of the wide orbit, consistent with the assertion that 
this triple system has a close to coplanar configuration. 



Key words: methods: data analysis - techniques: radial velocities - stars: 
binaries: spectroscopic - star: individual: HD 188753. 
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1 INTRODUCTION 



TODCOR (TwO-Dimensional CORrelation) is a two-dimensional correla tion technique to 
e radial velocities of the two components of spectroscopic binaries (jzucker fc Mazeh 
ODCOR was a natural extension of th e one-dimensional correlation technique 



1994). r 


rODCOR was 


(Griffin 


1967 




Simkin 





1974 



Tonry &: Davis 



19791 ) which turned out to be very effective 



for deriving radial velocities out of spectra with low signal-to-noise ratio (SNR). TODCOR 
calculates the correlation of an observed spectrum against a combination of two templates 
with different shifts, resulting in a two-dimensional correlation function, whose peak simul- 
taneously identifies the radial velocities of both primary and secondary components. Various 



studies have applied TODCOR 



lutions (e.g., 



Mazeh et al. 



1995 



successfully to many 



Metcalfe et al. 



of the m with quite faint companions (e.g.. 



binary spectra, vie 



1996 



Torres et al 



Mazeh et al 



1997 



1997 



ding double-li ned so- 



Zucker 



Prato et al 



20021 : 



2004T). some 



Mazeh et al 



2003|). 



Zucker et al. 



(119951 ) extended TODCOR to analyze spectra of triple-lined systems by cal- 
culating the correlation of an observed spectrum against a combination of three templates 
with three different shifts. The peak of the resulting three-dimensional correlation simulta- 
neously identifies the radial velocities of the primary, secondary and tertiary components. 
The new algorithm, TRICOR, was successfully applied to the spectra o f a few triple systems , 



yielding triple-lined solutions (e.g., 



Torres et al. 



1995 



JhaetaL 



2000). 



Torres et al. 



feooa) 



extended the algorithm to analyse spectra of quadruple systems. 

The importance of TRICOR arises from the prevalence of short-period b inaries accompa- 



nied by a distant comp anion, as is becoming evident from a few studi es (e.g.JMavor fc Mazeh 



19871 . and in particular 



Tokovinin fc Gorynya 



20071 ; see a review by 



Mazehl (12008 1). The ve- 



locities of the three components contribute to a better underst anding of the orbita l geometry 
of the triple system and the dynamics of the three stars (e.g., [Torres et al. 2006 ). 

Both TODCOR and TRICOR aim to analyze spectra composed of one order only. How- 
ever, many of the modern spectrographs produc e nowadays mult i-order spectra, some of 



Pepe et al. 



20001 ) . To optimize the use of 



which have as many as a few dozen orders (e.g., 
multi-order spectra for deriving accurate velocities and enabling us to detect faint compan- 
ions, we should in clude in the analysi s all the spectral information spread over the different 



orders. Therefore, 



Zucker et al. 



(120031 ) developed TODMOR, which can handle double-lined 



Multi-ORder spectra for the analysis of double-lined systems. Using TODMOR, 



Zucker et al. 
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(12003 



20041 ) could discover an evidence for a planet and a brown dwarf hidden in the spectra 



of the double-lined binary system HD 41004. 

The next step is to develop TRIMOR, a three-dimensional correlation function to handle 
multi-order spectra. This paper explains the basics of TRIMOR, which is conceptually easy 
to construct, once we have at hand TODMOR and TRICOR. In its second part, the paper 
applies the al gorithm to a set of spectra of HD 1 88753, which i s a well known triple system 



((Griffin 



19771 ) that attrac ted new interest lat e 



were obtained already by lEggenberger et al.l (120071 ) for a different study that focused on a 



Konacki 



20051 ). The spectra analyzed here 



careful analysis of the radial velocity of the brightest star in the system. Applying TRIMOR 
to those spectra enabled us to derive the velocities of the faintest companion, rendering the 
close pair of the system a double-lined binary. 

Section 2 presents the principles of TRIMOR, Section 3 applies TRIMOR to the already 
obtained spectra of HD 188753, while Section 4 derives the relative inclination of the system. 
Section 5 discusses briefly the results obtained for HD 188753 and future applications of 
TRIMOR. 



2 TRIMOR ALGORITHM 

2.1 The three-dimensional correlation 

Assume f(n) is an observed single-order spectrum of a triple system, where n is the wave- 
length. As f{n) is composed of three spectra, we wish to construct a model, m(n), by 
calculating a spectrum composed of gi, g 2 and g% — three templates for the primary, sec- 
ondary and tertiary components of the triple system, respectively. Given the templates, the 
model 

m(n) = gi(n - si) + ag 2 (n - s 2 ) + (3g 3 (n - s 3 ) (1) 

has five basic parameters: si, s 2 and S3, the Doppler shifts of the three templates, and a and 
/3, the flux ratios of g 2 and #3 relative to g±, respectively. The correlation — f(n) * m(n), 
between the observed spectrum f(n) and the model m(n), measures the goodness-of-fit of 
the model. Note that the correlation is not sensitive to any multiplication factor of the 
model function, and therefore only the relative flux of the three templates is relevant to our 
analysis. 
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TRICOR considers the correlation between the observed spectrum and the model as a 
function of all possible values of s\,s 2 and S3, generating a three-dimensional correlation 
function: 



-R (3) 0i> s 2 , s 3 ; a,(3) = f * [g x {n - s x ) + ag 2 (n - s 2 ) + (3g 3 {n - s 3 )] 



(2) 



where a and (3 are considered external parameters that have to be determined a priori. 
Alternatively, for each set of values of Si, s 2 and S3 we can find analytically th e values of a 



Zucker et al. 



(|1995h for details). 



and (3 that maximize the correlation (See Appendix A in 
The location of the maximum of R^ 3 \si, s 2 , S3) identifies simultaneously the radial velocities 
of the three components. To reduce the amount of computation we avoid direct calculation of 
the three-dim ensional corre l ation and instead calculate only six one-dimensional correlations 



as detailed in 



Zucker et al- 



ii^. 



To consider multi-order spectra we combine the correl ations o f the d ifferent orders ac- 
cording to the maximum-likelihood approach developed by IZuckerl (120031 ): 



M 



R 2 ( Sl ,s 2 ,s 3 ) = 1 - {IJ[1 -R-(si,s 2 ,s 3 )}} 



l/M 



(3) 



i=l 



where Ri is the three-dimensional correlation function of the i th order and M is the number 
of orders. 

When considering multi-order correlation functions, the values of a and (3 vary from 
order to order, and therefore, a and f3 are in fact vectors of values, {aj,i = 1,M} and 
{(3i,i = 1,M}. Finding numerically different a and (3 values for each order would greatly 
increase the number of fitted parameters. We, however, can reduce drastically the complexity 
of the calculation by taking advantage of the fact that for given two stars w ith known spe ctral 



Pickles 



19981 ). We 



types the flux ratio at a given wavelength is known approximately (e.g. 
therefore pre-calculated tables containing the fl ux ratio at ea ch order for many possible 
spectral-type pairs in our template library, using iPicklesl (119981 ) tables. In other words, the 
choice of templates determined the values of {a^i = 1,M} and {j3f,i = 1,M} in our 
analysis. 



2.2 Choosing the parameters of the correlation function 

To derive the most accurate velocities out of a given set of observed spectra we have to 
carefully choose the parameters of the model function. We first search for the best templates 
to fit the observed spectra. This step depends on the available library of templates that fit 
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the spectral range and resolution of the observed spectra. We first find the template for the 
primary, as this template has the largest contribution to the model function. We find the best 
template by running a one-dimensional correlation with each template candidate against all 
the observed spectra and choose the one that yields the highest correlation peak, averaged 
over all spectra. We then find the template for the secondary by running TODMOR with 
all template candidates over all observed spectra, and only at the third stage we choose the 
template for the tertiary. In many cases the contribution of the tertiary to the correlation 
function is minimal and we have to choose the template that leads to the best orbital solution 
of the tertiary, usually judged by the size of its residual RMS. 

As explained above, the choice of the three templates determines the flux ratio between 
the templates that we use. However, to maximize the correlation we allow a multiplicative 
factors f a and fp of the order of unity which we apply to the pre-determined flux ratios of all 
orders. The factors are determined numerically to maximize the average of the correlation 
peak. 

In some cases the actual stars in the system have rotational broadening larger than 
their templates. To correct for that we try to insert numerically some additional rotational 
broadening to the three templates. This is done by convolvin g the template spectr a with 
a broadening function, based on the Fourier transformation of ISmith fc Grayl (119761 ). Here 
again we first determine the rotational broadening of the primary template, then the sec- 
ondary and at last the tertiary. The values of the rotation of the three templates are chosen 
to maximize the averaged correlation peak. 

In finding the best parameter values of the model function, we always choose the ones 
that maximize the correlation peak averaged over all observed spectra. Because the values 
of all these parameters are based on the whole set of observed spectra we consider these 
parameters as global ones. In some cases, however, we prefer to perform the averaging only 
over the high-SNR spectra, as the low-SNR spectra add noise to the process. 

To conclude, we fit eight global parameters to the whole set of observed spectra: the 
best three spectral templates, two multiplicative factors to be applied to the flux tables, 
and three rotational broadening velocities. After finding these eight global parameters, we 
consider the three-dimensional correlation for each observed spectrum, find its maximum 
and derive the three velocities of the three stars for each spectrum. 

The uncertainties of the three velocities were derived as explained by IZuckerl (120031 ). by 



deriving the inverted Hessian of the likelihood of the velocities. However, these errors were 
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proven to be overestimated in quite a few cases (e.g., lEggenberger et al.l 120071 ). and we had 
to rescale them by using the RMS of the residuals of the orbital solution. In this case, we 
kept only the relative sense of the derived error estimate, which probably depended mainly 
on the SNR of the obtained spectra. 



3 APPLICATION OF TRIMOR TO THE CASE OF HD 188753 



3.1 The system HD 188753 



HD 18 8753 (H I 



1999 . 



\. iHoughl 



980 01) is a bright triple system, 7.4 V-mag, at a distance of 46 pc (jSoderhjelm 



Hough! (118991 ) was the first to discover that the system was a visual binary, with a 



separation between the two components, A and B, of ~ 0.3 arc-sec. The AB pair is char- 



acterized by a per iod of 25 .7 years, a sem i-major axis of 



(jSoderhjelm 



1999, see also 



Konacki 



2.3 A U and an eccentricity of 0.5 



20051 ) . Later iGriffinl (119771 ) discovered that the system 



contained a spectroscopic binary with an orbital period of 155 days and an eccentricity of 
0.26. 

Recently, HD 188753 attracted interest when iKonackil (120051 ) reported a radial- velocity 



modulation of the primary — A, with a period of 3.35 days and an amplitude that implied 
an unseen companion with a minimum mass of 1.14 Jupiter mass. The detection of a close-in 



Jupiter-like planet in the trip 
planet formation theory (e.g., 



e system HD 188753 posed a p roblem for the widely accepted 



Pfah 



2005 



Jang-Condell 



close-in planet is now a matter of debate, after 



2007 ). However, the existence of the 



Eggenberger et al.l (120071 ) obtained their own 



ELODIE spectra, the analysis of which showed no evidence to support the conjecture of a 
pla net in the system. 



Eggenberger et al.l (120071 ) used TODMOR to derive the velocit ies of A, the bright distant 



Eggenberger et al. 



(|2007h 



star, and Ba, the brighter star of the close pair. The main goal of 
was to study the radial velocity of component A, and to confirm or refute the planet con- 
jecture. Nevertheless, their spectra included information on the close pair, which before this 



work remained a single-lined spectroscopic binary. To demonstrate the 



capab ility of TRI- 



MOR we re-analysed here the spectra obtained by lEggenberger et al.l (120071 ) and derived 
the velocities of A, Ba and Bb, the latter being the faintest component of the triple system, 
rendering the close binary a double-lined spectroscopic system. 
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3.2 The TRIMOR radial velocities 



The data we ana 



Eggenberger et al. 



yzed h ere comprised the ELODIE spectra of HD 188753, obtained by 
(120071 ) between July 2005 and March 2006. ELODIE is a fiber-fed spec- 
trograph at the Obs ervatoire de Haute-Provence (France) with a fiber of 2" in diameter 
( IBaranne et al.lll996l ). With a resolution of A/ A A = 42 000 it covers the wavelength range 
3850-6800 A with 67 orders. The analysis used 32 ELODIE orders within the spectral range 
4810-6800 A, after having excluded orders that are heavily polluted by telluric lines. We 
also excluded the blue orders because in these orders we expected the secondary and ter- 
tiary signals to be too weak relative to the signal of the hotter and bluer primary. 

Our library of templates included high-SNR ELODIE and CORALIE spectra, with spec- 
tral types that run from F8 to M6. To choose the best configuration of templates and ro- 
tational broadening we proceeded as detailed above. In our case, the tertiary template was 
so faint that its impact on the correlation peaks was minimal. Therefore, we carried out 
the analysis with the different possible tertiary templates through the orbital solution, and 
chose the one that minimized the residual RMS of the tertiary orbit. 

We finally chose a G6 dwarf (HD 224752) for A, a K0V dwarf (HD 225208) for Ba, and 
a K7 dwarf (GL 338B) for Bb. As we do not have a table with flux ratios for every pair of 
spectral types, we used for a our G5K0 table, for which the flux ratio at 6800 A was 0.53, 
and a multiplication factor, f a , of 1.12. For (3 we used the G8M0 table, for which the flux 
ratio at 6800 A was 0.06, and a multiplication factor, fp, of 0.91. In our model the fluxes 
of the primary, A, and the secondary, Ba, are comparable, while the tertiary, Bb, is much 
fainter than the other two. The additional rotational broadening velocities chosen for the 
templates used were 0.0 kms -1 , 2.0 kms" 1 and 2.5 kms -1 for A, Ba and Bb, respectively. 



As in 



Eggenberger et al.l (120071 ). we also noticed that when two of the three stars had 



similar velocities, TRIMOR sometimes picked the wrong peak of t he correlation, yield i ng er 



roneous velocities. Three additional spectra that were analysed by lEggenberger et al.l (120071 ) 
yielded ambiguous results. We preferred not to include these velocities in our analysis. Thus 
we ended up using only 35 spectra, with a typical SNR of 55 per pixel at 5500 A. 

Our radial velocities and their error estimate for the three components of HD 188753 are 
given in Table [U The original error estimates for the velocities of Ba and Bb turned out 
to be larger than the actual residuals of the orbital s olution. This effect appea red in a few 
other stellar systems that we have analysed (see also lEggenberger et al.l (120071 1 for similar 
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Table 1. The derived radial velocities for HD 188753 



JD-2450000 Primary [kms x ] Secondary [kms 1 ] Tertiary [kms 1 ] 
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effect with TODMOR). We therefore adjusted the error estimates by multiplying the derived 
errors of Ba by 0.34 and those of Bb by 0.26, so the normalized x 2 value of the solution is 
unity for each of the three stars. 



3.3 The orbital solution 

To derive the orbital solution we solved for double-lined solution of Ba and Bb, simul- 
taneously with a linear drift of their center-of-mass velocity. This added one more linear 
parameter to the function fitted to the radial velocity measurements. For the A velocities, 
we solved independently for a linear drift only. 

The orbital solution based on the TRIMOR velocities is given in Table [2] and plotted in 
Figure HJ The linear drift for A and B could barely be detected, because of the relatively 
short time span of the observations. All our derived values are consistent with the results 
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HD 188753 




3600 3650 3700 3750 3800 3850 3900 
JD-2450000 (days) 



Figure 1. Orbital solution for HD 188753. The circles refer to A, triangles to Ba and squares to Bb. Error bars in A are smaller 
than the size of the data points. The solid straight line represents the 25.7-year orbital motion of A and the dashed line that 
of B. 



Table 2. Orbital parameters for the 154-day spectroscopic orbit of HD 188753 B. The long-period orbital motion of the AB 
pair was taken into account by including a linear drift. 



Parameter 


Units 


Value 


P 


(days) 


154.448 ± 0.015 


r 


(JD-2400000) 


53405.007 ± 0.032 


e 




0.175 ±0.001 


7 


(kms" 1 ) 


-21.623 ± 0.014 


Li-' 


(dug) 


134.87 ±0.02 




(kms" 1 ) 


13.479 ±0.017 


K 2 


(kms" 1 ) 


17.63 ±0.09 


Linear drift A 


(km s — 1 yr — - 1 ) 


0.619 ±0.056 


Linear drift B 


(km s" 1 yr — 1 ) 


0.346 ± 0.039 


mcas 




35 


rms (A) 


(ms- 1 ) 


98.3 


rms (Ba) 


(m s -1 ) 


72.1 


rms (Bb) 


(ms — - 1 ) 


455.6 



of 



Eggenberger et al. 



Eggenberger et al. 



(|2007t) . Our error estimates of these values are smaller than those of 
(120071 ). indicating another possible advantage of TRIMOR. 
The residuals of the three sets of velocities are plotted in Figure El While the residual 
RMS of A and Ba is at the level of 0.1 km s _1 , the RMS of the Bb residuals is of the order 
of 0.5 km s _1 . This is consistent with the relative brightness of the three stars. 
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Figure 2. The residuals of the orbital solution in Fig. [T] The RMS for A is 98 ms , for Ba is 71 ms 1 and for Bb is 0.5 
km s — 1 . 



The new double-lined solution of the B system yields directly its mass ratio: 
m Bb 



0.768 ± 0.004 . 



m Ba 



The mass ratio of the wide orbit is: 

w,b 

q AB = — = 1.79 ±0.26 . 

m A 



(4) 



(5) 



The latter has a large error, becaus e of th e short time span of the observations, but still is in 



agreement with 
way. 



Eggenberger et al. 



(120071 ). who derived the linear drift of A with a different 



3.4 Search for evidence of a planet 

To search for an evidence of the planet suggested by 



Konackil (120051 ). we performed a search 



Eggenberger et al 



for pe riodicity in our derived radial velocities of HD 188753 A, as did 
(120071 ). In fact, the residuals of A were slightly larger than those of B, and therefore could 
hide additional real variabilit y. To search for any periodicity we derived a double-harmonic 
periodogram of the residuals (jShporer et al.l 120071 ). by fitting for each trial period a double- 
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8 - 




Frequency (days 1 ) 



Figure 3. Double-harmonic periodogram of the residuals around the linear trend for HD 188753 A. For this data set, the 1% 
false alarm probability corresponds to a power of 8.58 and is represented by the top of the box. The dashed line denotes the 
frequency of the planetary signal reported by Konacki. 



harmonic function. The periodogram value for each trial period was taken as the power (sum 
of the squared coefficients of the two harmonics) divided by the x 2 goodness-of-fit parameter 



of the fitted 



(.e.g. 



Scargk 



ight c urve. We preferred this analysis to the frequently used Lomb-Scargle test 



19821 ). as the radial velocity of the presumed additional planet might not have a 



perfect circular orbit, and therefore might have a periodic modulat i on wi th a few harmonics. 
Our method is a further generalization of Izechmeister &: Kursteil ( 2009 ) method. 

To estimate the statistical significance of a peak with an height S we analysed 10,000 
randomized set of residuals and found the fraction of periodograms that had a peak higher 
than S. The randomized residuals were derived by taking the actual residuals, wit h the actual 
timin gs randomly permutated. Assuming a negligble red noise in the data (e.g 



Pont et al 



20061 1 . this test gives the statistical significance of any detection. 



The computed periodogram of the residuals around the linear trend for component A 
are plotted in Figure [3j Our velocities for HD 188753 A show no sign of a short-period signal 
with the period of ~ 3 days, although some variability with a period of about 20 days might 
be present in the data. However, the periodagram peak might reflect the window function of 
the data, and in any case is not significant enough to draw any conclusion without further 
observations. 
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4 THE RELATIVE INCLINATION OF THE TWO ORBITS 

HD 188753 is a special case, for which accurate astrometric, photometric and spectroscopic 
information is available. We wish to use this information and estimate the degree of alignment 
between the or bital planes of the wide (astrometric) and the close (spectroscopic) orbits. 
Konackil (120051 ) shortly addressed this issue, but he did not have yet a direct measurement 
of t he mass rat i o of th e close pair. 



Soderhjelml (119991 ) estimated the inclination of the wide orbit at about 34°. We need 



now to derive the inclination of the spectroscopic binary so we can estimate the relative 
inclination. 

For any double-lined spectroscopic binary, the orbital inclination can be written as: 



sin i 



4.7 x 10~ 3 VT^ 



K x + K 2 



P 



1/3 



'M t + M 2 



-1/3 



(6) 



lkms~V V lda y/ \ M Q J 
In this formula, only the total mass of the spectroscopic binary is not known from the radial- 
velocity data. Deriving the total mass of the binary from information not included in the 
orbital solution leads to an estimate of the binary inclination. 

We will use here three approaches to estimate the B-component mass. As in many SB2s, 
we can derive the total mass, or at least estimate it, through analyzing the spectra of the two 
stars. In our case we can take advantage of two additional known features of HD 188753. 
First, the binary is a component of a wide orbit with known relative astrometric orbit. 
Therefore, the total mass of the triple system is known. To derive the mass of the close pair we 



only need the mass ratio of the wide orbit, w hich we can obtain 



provided their time span i s long enough (e.g . 



of the system is known (jSoderhjelm 



Eggenberger et al 



rom t he spectroscopic data, 
2007j ) . Second, the parallax 



19991 ). because A is bright enough to be measured 



by Hipparcos. Using the parallax, we can derive the expected magnitude of HD188753Ba 
and HD188753Bb for each assumed mass of the Ba component, using the mass ratio, qs, 
obtained from the radial-velocity solution. We can then compare the expected magnitude 
of the combine d light of Ba an d Bb with the observed J, H and K s 2MASS magnitudes of 
B obtained by iKonackil (120051 ). In this way we can find the actual mass of Ba, the single 
remaining independent parameter in this exercise. 

Table [3] presents the results of the three approaches. For the spectral analysis approach 
we utilized the spec tral types we used in o ur TRIMOR analysis (KOV and K7V) and a 
standard calibration (IHabets fc Heintzdll98ll ) to convert th em to masses. For the astrometric 



approach we used the total mass of the system, 2.73M (jSoderhjelml Il999l ) . and the mass 
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Figure 4. For each assumed mass of the Ba component, the plot shows the derived K B magnitude (in the 2MASS K B band), 
assuming the known mass ratio qg, and the known distance modulus of 3.31. The dashed line represents the observed K B 
magnitude. 



Table 3. The total mass of the close binary and its inclination, based on three approaches 



approach 


M Ba + M Bb 


sin i 


% 




[Mq] 






spectral types obtained by TRIMOR 


1.40 ±0.15 


0.690 ± 0.026 


43?6±2?1 


astromctry of the visual pair 


1.78 ±0.31 


0.637 ±0.038 


39?6±2?8 


2MASS photometry 


1.634 ±0.019 


0.6557 ± 0.0079 


40? 97 ±0? 60 



ratio derived by lEggenberger et al.l (120071 ). For the p hotometric approach w e used a parallax 



of 7r = 21.9 ±0.6 mas, the pho tometric calibration of 



Bessell fc Brett 



(119881) and the 2MASS 



correction of ICarpenterl (1200 ll ). Figure H] demonstrates this procedure for the K s band. Error 
estimates are not readily available for the photometric calibrations, but we can use the 
variance among the results from the three photometric bands as a proxy to the error. 

The results of the three approaches show a remarkable consistency, indicating that the 
inclination of the close pair is around 41°. It is not easy to estimate the uncertainty of the 
derived inclination, as the biases of the approaches are not well known. Arbitrarily, we assign 
an error of 3° to our estimate of the binary inclination. 

The close-pair inclination found here, 41° ± 3°, is close to 34°, the value given by 
(jSoderhjefml 119991 ) for the wide orbit inclination. This means that the minimal relative 
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inclination (MRI) between the two orbits is around 7°. ISoderhjelml (119991 ) values are not 
accompanied by error estimates, and therefore it is difficult to assign an uncertainty to the 
MRI value of the triple system. However, if we arbitrarily assign an error of 3° to the wide 
orbit inclination, we can deduce an error of ~ 5° to our MRI. Our findings are therefore 
consistent with the MRI being close to zero, and therefore the system might have coplanar 
configuration. 



5 DISCUSSION 

In this paper we laid out the principles of TRIMOR, an algorithm to analyse multi-order 
spectra of triple systems and deduce the radial velocities of the three components. A few 
improvements of the algorithm could still be applied, one of which is the derivation of the 
flux ratios of each order. Instead of relying on pre-determined ratios, based on the template 
spectral types, one could derive the flux ratios from the observed spectra, as done in TOD- 
COR and TRICOR. The best flux ratios could be derived by maximizing the correlation, 
using all observed spectra and therefore can, with high enough SNR, yield better results. 
Another improvement, also associated with the flux ratio, is to apply a ratio that is varying 
along each order. That is, instead of the model of Equation (1) for each order, we can use a 
modified model: 

R {3) (s u s 2 , s 3 ; a(n),P(n)) = f*\gi(n- si) + a(n)g 2 (n - s 2 ) + P(n)g s (n - s 3 )] . (7) 

The functions a(n) and f3(n) could be approximated by low-degree polynomials. In fact, 
such an improvement can be applied to TODCOR, TRICOR and TODMOR as well. 

We have applied TRIMOR to the obtained spectra of HD 188753. Without investing 
any additional observational resources we derived the velocities of the faint component of 
the triple system and turned the close pair into double-lined spectroscopic binary. The new 
radial velocities led to an estimate of the close pair inclination, which depends on the total 
mass of the close pair. We estimated the total mass of the close pair by three approaches, 
one that relied on the spectral type of the two stars of the close pair derived here, one that 
used the astrometry of the wide orbit and the mass ratio of the wide orbit, independent of 
the present analysis, and finally another approach that depended on the mass ratio of the 
close pair derived from the new velocities of the present analysis and known photometry 
of the close pair. All three approaches yielded similar close-pair inclinations that are close 
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to 41°. The consistency of the three estimates of the inclination demonstrated the power of 
TRIMOR algorithm and established more confidence in the present result. 

The derived close-pair inclination suggests that the relative inclination of HD 188753 is 
close to zero and the motion is close to being coplanar. This result adds HD 188 753 to the 



Tokovinin 



2004 



Tokovinin et al. 



2006 




Tokovinin 


2008). 


r 


Shkolnik et al. 


200s| 



2008). Applying TRIMOR to more close triple or 
for which we have at hand the astrometric 



quadruple systems (e.g., 
orbit of the wide pair from Hipparcos, can enlarge substantially the set of triples with 
known relative inclinations. This certainly will be true after the new GAIA satellite all-sky 
astrometry will be available in the near future (jPerryman et al.l 120081 ) . 

The distribution of the relative inclinations in triple systems can shed s ome light on 



Sterzik fc Tokovinin 



formation models of binary and triple systems (e.g.. 
be relevant to the tidal evolution of close multiple system s (e.g., 
multiple system formation m odels of fragmen t ation (e.g., 
later stage of accretion (e.g., 



Bate et al 



2002 



Mazeh 



Machida 



Bonnell fc Bat 



2002) and can also 



20081 ). Binary and 



20081 ). which include a 



20061 ). predict alignment of 



the rotation stellar axes with the orbital angular vector in binaries, as well as coplanar 
motion in triple systems. On the other hand, m odels assuming that binaries and multiple 



Clarke 



1996 



Delgado-Donate et al J 12004 ) 



systems are formed by N-body interaction (e.g., 
predict that the stellar rotation might not be aligned with the orbital motion, and triple 
systems could have high relative inclination. 

The relative inclination of triple systems might also have implication for long term evo- 
lution of systems for which tidal evolution of t he clo se pair can circularize the binary orbit 



during its stellar life time. 



Mazeh fc Shahaml (119791 ) found that the third star can induce 



almost periodic eccentricit y modulation into the binary motion, even if the binary starts 
with zero ecce ntricity. (See iKozai (119621 ). who discovered this effect for asteroids with initial 



eccentricities.) iMazeh &: Shahaml pointed out that the eccentricity modulation implies that 



binaries with close enough distant companions never get into a stable circular orbit, which 
means that the binary tidal dissipation continues to be active forever. The injected binary 
eccentricity invokes frictional forces inside the two stars that dissipate rotational energy and 
transfer angular momentum between the stellar rotation and the binary orbital motion. This 
causes the semi-major axis of the binary to shrink on the circularization time scale. 

The possible impact of the third star on the sh rinking of the close pair se paration depends 
on the amplitude of the eccentricity modulation. IMazeh fc Shahaml (119791 ) showed that this 
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amplitude does not depend on the third dis tant star separation, but strong l y dep e nd on the 



relative inclination between t 



re two orbits. 



(120061) 



Borkovits et al 



found by 



Mazeh fc Shahaml ( 



Eggleton fc Kiseleva-Eggletonl (j200ll ). Eggleton 



(2004 ), an d lFabrycky fc Tremaind (120071 ) further studied the effect 



19791 ). recently termed as the KCTF (Kozai Cycle with Tidal 



Friction) effect, and showed that the tidal effect that leads to the shrinking of the close- 
pair separation occurs when t he relative angle betwe en the two orbital planes of motion 



is larger than ~ 40°. In fact, 



Tokovinin et al. 



20061 ) suggested that all or at least most 



of the very close binaries are found wi th a distant compa nion, sug gesting that the y were 



formed through the effect suggested by 



Mazeh fc Shaham (see, also 



Tokovinin 



2008 



for an 



updated statistical analysis). It is therefore important to get the true distribution of relative 
inclinations of triple systems, with close pairs in particular. Cases like HD 188753, with small 
relative angle, do not support the shrinking scenario of close binaries in triple systems, as the 
eccentricity modulation expected in such cases is small. As we have shown here, TRIMOR 
can be very effective in obtaining the relevant information. 

TRIMOR can also be effective in searching spectroscopic evidence for unknown faint third 
companions in spectra of known double-lined binaries. The accumulated data base of multi- 
order spectra of binaries can be use d for such a proj e ct. A somewhat similar approach, with 



a different algorithm, was taken by 



D'Angelo et al.l (120061 ) . who searched for spectroscopic 



trace of a faint tertiary in a large set of observed spectra of contact binaries. In about a third 
of the systems they found an evidence for a tertiary. We note that their spectra are composed 
of one order centered on 5200 A, a wavelength for which K and M stars, for example, are 
quite faint. In our case, where we had the luxury of using CORALIE spectra whose orders 
get up to 6800 A, we could detect and measure radial velocities of a late K-type tertiary star. 
The search for faint companions could be even more e fficient with spectra that go further 
towards the red wavelegths (e.g.. iFabricant et al.ll2008l ). 



Finally, TRIMOR can help find extra-solar planets in triple systems, in t 



re same way 



TODMOR was instrumental in discovering the planet in the HD 41004 system (IZucker et al. 



2003 



200J). 
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